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On , Abstract: We study how topological defects manifest themselves in the equal-time two- 

point field correlator. We consider a scalar field with Z2 symmetry in 1, 2 and 3 spatial 
^D ' dimensions, allowing for kinks, domain lines and domain walls, respectively. Using nu- 

10 ' merical lattice simulations, we find that in any number of dimensions, the correlator in 

momentum space is to a very good approximation the product of two factors, one describ- 
ing the spatial distribution of the defects and the other describing the defect shape. When 
the defects are produced by the Kibble mechanism, the former has a universal form as a 
k> I function of k/n, which we determine numerically. This signature makes it possible to de- 

H ' termine the kink density from the field correlator without having to resort to the Gaussian 

approximation. This is essential when studying field dynamics with methods relying only 
on correlators (Schwinger-Dyson, 2PI). 
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1. Introduction 

Spontaneously broken global and local symmetries play an important role in many physi- 
cal systems, from high energy physics to condensed matter. In many cases, the symmetry 
breaking pattern is topologically non-trivial, giving rise to a possibility of topological de- 
fects, non-linear objects which are stable because of their topology. Examples of these 
include domain walls in ferromagnets, vortices in superfluids and superconductors, and 
magnetic monopoles and cosmic strings in high energy physics and cosmology. 

In general, spontaneously broken symmetries are restored at sufficiently high temper- 
atures. As the system cools down, it undergoes a phase transition into the broken phase. 
When this happens, the system locally picks out one of the possible vacuum states. This 
choice is dynamical and random, and would generally have different outcome in different 
patches separated by more than one correlation length. This phenomenon, known as the 
Kibble mechanism, results in the creation of topological defects |j|, ^, ^. It has been stud- 
ied in detail in many experimental setups Q, |5|, |^, 0, ^, ^, 10, 11, 12, 13 1, and it may have 



also taken place in the early universe, producing domain walls, cosmic strings or magnetic 
monopoles. 

The number of defects created in such a transition depends on the dynamics of the 
model, the nature of the symmetry and the cooling rate, which in the cosmological context 
is determined by the Hubble rate H. This dependence can be estimated by considering the 
critical behaviour the theory and expressed in terms of its critical exponents 0] , and these 
estimates have been shown to be accurate in some condensed matter experiments ||9|, |lO| 
and lattice field theory simulations |14, [I^, ^, ^]. In these simulations, the classical 



equations of motion are solved numerically on a spatial lattice, and the number of defects 
was determined by identifying them in the final field configuration and counting them. This 



approach is fully non-perturbative and can be used also in gauge field theories |18, 19], but 
it cannot incorporate quantum mechanical efi'ects and even in a classical theory, any noise 
such as thermal fluctuations can make the defects hard to distinguish, and the counting ill 
defined. 

Because of these limitations, it would be advantageous to study defect formation from 
first principles using other techniques that would give a more detailed understanding of 
the process and would be valid also in quantum theory. There have been many attempts 
to do this using various techniques ranging from a linear approximation |20, ^, 22, ^] to 



Schwinger-Dyson equations, Hartree approximation [24, 14] and the 2PI formalism ]25]. 

A common feature in these alternative approaches is that one does not have access 
to individual field configurations, and therefore direct counting of defects is not possible. 
Instead, the dynamical variables are correlation functions. This raises the question of how 
to determine the number of defects from the field correlator. If one assumes that the field 
ensemble is Gaussian and identifies zeros of the field with defects, one can derive a simple 



expression for their density in terms of derivatives of the two-point correlator ]^, 27, ^, ^ . 
However, because defects are non-linear objects, their presence makes the field ensemble 
non-Gaussian, and therefore this assumption is not justified. Indeed, numerical simulations 



have shown that this approach does not work ]29]. 

In this paper, we show how kinks, domain lines and walls in 1, 2 and 3 spatial dimen- 
sions respectively, manifest themselves in the equal time two-point correlator in momentum 
space at a non-linear level, without relying on the Gaussian assumption. This is an ex- 
tension of previous work ]|2^, in which we showed that the 2PI formalism is unable to 
reproduce the classical signatures of global defects, at least in a 1/A^ expansion to next-to- 
leading order. 

We will first set up a simple model describing kinks, domain walls and domain lines 
(section |2|), and then derive the signatures in the two-point function which we will be 
looking for (section ^). We then perform sample lattice simulations in both 1, 2 and 3 
spatial dimensions (section 0), and demonstrate the signatures in practice, calibrating the 
net defect density the a benchmark density (section S) . We conclude in section pi 

2. Setup, model and defects 

We consider the action of a real scalar field (j) in D spatial dimensions 

s= I d'^xdt {^d^^d^'(p-v{(i))\ , (2.1) 

with a potential V^(</>), so that the theory has Z2 symmetry {(f) o —(/)), 

n'^) = -^MV + ^0^ (2.2) 

The parameters and fields have energy dimensions [0] = {D — l)/2, [^] = 1, [A] = 3 — D. 
The equation of motion reads, 

d'i(t>{x,t)+Tdt(t>{x,t) - dl(t>[x,t) - /i^(/'(x,t) + -0^(x,t) = 0. (2.3) 



and we have added a small damping term Tdfcl) ([r] = !)• This will drive the system from 
a high temperature initial state to a cold, near- vacuum final state (see also section |^. 

The model has two degenerate vacua at ipo = if = ±y^6/X ^, and there are topological 
defects, kinks {D = 1), domain lines {D = 2), domain walls {D = 3), which interpolate 
between them. The classical kink solution is 

X 

(/'kink (a;) = f tanh-, (2.4) 

where d = v2/^ is the kink thickness. Domain lines and domain walls are extensions in 1 
and 2 dimensions of kinks, possibly with curvature on a length scale which we will assume 
to be much larger than d. 

In the following, we consider a situation in which a network of defects has been formed 
in a symmetry breaking phase transition by the Kibble mechanism. In this case, the defect 
network has one characteristic length scale, which is determined by the field correlation 
length at the time of the transition. It corresponds to the typical distance between defects 
or, equivalently, the number density n of defects. 

3. The correlator ansatz in 1, 2 and 3 dimensions 

Our basic observable is the two-point correlator 

G{r = \x-y\) = {(t>{x)(t>{y)), (3.1) 

assumed here to be homogeneous and isotropic, and therefore a function of the relative 
position r only. The only two scales in the system are the kink thickness d (or /u), and 
the defect density n^. In general, we would expect the position and momentum space 
correlators to be some general functions G{r,d,n) and G{k,d,n). In [^], we argued that 
in ID, in the limit of (i = 0, and if kinks are distributed at random with average density 
n, the correlator should have the form 

G(.,„) = .^e--, G(t,„) = |^^^|^. (3.2) 

If the kinks all have the same smooth profile 4'\dn\s.{x) ^ this modifies the correlator by a 
multiplicative factor, 

where (/>kink(^) is the Fourier transform of the kink profile (ET 



2w lirkd , ,^ 

'^kink(A:) = ^-^TT-ry (3-4) 



k sinh -^'Kkd 



^Related in a complicated way to the damping rate F. 



Therefore we should find 

G(k,d,n) = ,, , ,A ^TH r . 3.5 

"■ ' ' ' 2nl + {k/2nf ysmh{lTTkd) J ^ ' 

Since d and v are known^, this expression allows for a one-parameter fit in n. More 
generally, the field correlator factorises as 

G{k,d,n) = —GcorT{k/n)G],ink{kd), (3.6) 

n 

where Gcorrik/n) describes the spatial distribution of the kinks and Gkink(^^) the kink 
shape. In particular, when the kinks are formed by the Kibble mechanism, the only relevant 
scale is their density, and therefore Gcorr is a function of the dimensionless combination 
k/n only. 

We will generalise this procedure and instead make the ansatz that the field correlator 
factorises in the same way in higher dimensions. 



V 



2 



G{k,d,n) = -^G^^„{k/n)GM,,]^{kd). (3.7) 

Several comments are in order at this point. Firstly, the d and n dependence are separated, 
and apart from the trivial scaling n~ coming from the measure in the Fourier transform, 
the correlator only depends on the dimensionless quantities kd and k/n. Also, we can 
obtain the kink correlation function, as 

^ ,, , , n^ Gik, d,n) , , 

The Kibble mechanism predicts that this is a universal function of the ratio k/n only, and 
therefore if we measure this for different parameters and different cooling rates, the results 
should coincide. We will use numerical simulations to demonstrate that this is indeed the 
case, and that to a very good approximation, the universal function (in coordinate space) 
is of the form (for some a^, 6j) 

Gcorr(nr) = (aie-«2("'-)' + b,e-'^^^^ , (3.9) 

which in D-dimensional momentum space becomes (for some other ai, I3i) 

G£„ (k/n) = ( ai e-°^ ('^Z")' + ^' m^iva ) • (3- 10) 

V [l + /32(fc/n)2](^+i)/2y 

We use our numerical results to determine the parameters a^ and /3j. After this, the form of 
the correlator is fixed, and the only remaining free parameter is n which sets the scale. The 
defect density n can therefore be determined directly from the correlator using a simple 
one-parameter fit. 



^Strictly speaking, these are the vacuum values, but we will see that they are sufficient. We will comment 
on finite temperature effects below. 



4. Numerical procedure 

To carry out the numerical calculation, we discretised the equation of motion ( |2.3| ) on a 
lattice of spacing a and using a leapfrog algorithm with time step a6t for the time derivative. 
Input parameters are lattice size n^, St, and the dimensionless o/i, aT, a^~^X. 

We choose the initial conditions at time f = to mimic the quantum vacuum state 
corresponding to the potential 

yini(0) = ^fl''^^. (4.1) 

Because this is a free theory, the equal-time quantum two-point functions of the field (j) 
and its canonical momentum vr = dt(j) are simply 

{Hk)m) = (2vr)5(A: + 'Z)^, (vr(A:)^(g)) = (27r)<5(A; + g)^, (<A(fc)7r(g)) = 0, (4.2) 

where u^ = \/k'^ + M^- Our initial condition is realised by a Gaussian ensemble of field 
configurations which has these same two-point functions |3^. We initialise all modes, and 
although this corresponds to a divergent energy in the continuum limit, at finite lattice 
spacing it is just a specific choice. In our case, it is not crucial how well these initial 
conditions reproduce the actual initial quantum state, since we are only interested in the 
classical dynamics. For our purpose, it is sufficient that it corresponds to a state with 
unbroken symmetry^. 

The classical equation of motion ( |2.3D allows us to rescale the coupling A to unity, sug- 
gesting that only the dimensionless ratio A//i^ plays a role. However, the initial conditions 
( |4.2| ) remove this freedom. 

We solve the time evolution of the system using the classical equation of motion ( p.3| ) 
for a large number of initial configurations that are drawn from the distribution specified 
by Eq. ( |4.2D . Because of the inclusion of a damping term F, the system will lose energy and 
cool down to a state close to the vacuum. However, as we are going through a symmetry 
breaking transition, topological defects will be created. At the end of the evolution, these 
will be "frozen" in, possibly with very small residual thermal fluctuation superposed. The 
correlator at this final time is our observable, from which we will extract the number of 
defects. 

5. Calibration and results 

To determine the parameters a^, /3j in ( |3.10D from the numerical data, we need to have an 
independent measurement of the defect density n. We do this by direct counting of zeros. 
We count, configuration by configuration, all lattice points A'^points where the field value is 



^ Other possible choices could be a Bose-Einstein distribution (finite UV energy), or a classical thermal 
distribution (divergent UV energy). We found that it made very little difference to the final number of 
defects, as long as the initial energy density was high enough. See for instance [ p2| , B3, p4j for studies of 
different initial condition prescriptions in the context of symmetry breaking and thermalisation. 
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Figure 1: The "counting" defect density n for different F in 13 = 1,2,3 dimensions (circles). 



Superposed, the density n determined from a fit of the field correlator to our ansatz (3.10) 



in some range [— ^, 6\ around </> = 0. We divide by the width of the kink profile in that 
range, 



tanh[Ax/(i]^ < 6*, d= V2/{an). 



giving for instance 



e = 0.1, -^ Ax{e) = ±0.327ri, 9 = 0.05, -^ Ax{9) = ±0.227(i, 
0.025, -^ Ax{9) = ±0.159(1, 9 = 0.0125 -^ Ax{9) = ±0.112(i. 



(5.1) 



(5.2) 
(5.3) 



The number of kinks/length of wall is then the number of lattice points divided by this 
length, and the density (in lattice units) is 



ni9) 



N, 



points 



2Ax{9)n^N, 



configs 



(5.4) 



This we consider to be the benchmark defect density, at least for late times. We checked 
that except at very early time, when the walls are not well separated, the kink/wall density 
is practically independent of the threshold 9. This conclusion holds for D = 1,2,3. The 
counting density as a function of damping rate is shown in Fig. ||. 

For a given F, we calculate this "counting" density, and use it to rescale the correlator 
as in Eq. (p.8|). If the resulting curves coincide, we will have shown that 1) the original 
correlator is separable in kd and k/n, 2) the dependence on n is through nr only 3) n is 



in fact the "counting" n (5.4), up to an overall calibration and 4) the rescaled curve is 
Gcorr) which One can then attempt to approximate. Note that since the density ranges over 
orders of magnitude (Fig. ||), the curves certainly do not coincide without rescaling. 

Fig. g shows the rescaled curves for various F for simulations in D = 1 spatial dimen- 
sions. We use a lattice of 2^^ = 131072 points, afi = 0.1, and a^X = 0.6. We have f^ = 0.1, 
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Figure 2: The rescaled propagators (3 
in D = 1 for different F. 



• Data 

— Gaussian component 

— Exponential component 

— Fit 




Figure 3: One rescaied propagator (F/yit = 
0.1) witlr Irigti statistics. Superposed, tfie fit- 
ted form ( |3.9| ), as well as the Gaussian and 
exponential components separately. 



d = 1/2/ a// ~ 14.14. We consider a range of damping rates ^ = 0.05 — 3.2. The curves 
agree very well indeed. In Fig. ^ we show a single rescaled curve with improved statistics 
(16 times as many configurations) and the overlaid fit of the anticipated form (3.9). The 
fit is quite convincing, although at intermediate scales, there is an oscillating feature which 
we choose not to account for in our description. The fit parameters Oj, /3j are given in 
Table |l[ 



D 



ai 



0L2 



/?: 



/32 



1 0.384 ± 0.001 

2 0.725 ± 0.008 

3 1.73 ±0.01 



0.0356 ± 0.0001 
0.1039 ±0.0006 
0.1654 ±0.0004 



0.180 ±0.001 
0.274 ±0.010 
0.406 ± 0.013 



0.0447 ± 0.0004 
0.097 ±0.003 
0.110 ±0.003 



Table 1: The fit parameters in (3.1C) determined from the numerical data. 
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Figure 4: The rescaled propagators (| 
£) = 2 for different F. 



■ Data 

— Gaussian component 

— Exponential component 

— Fit 




Figure 5: One rescaled propagator (F//i = 
0.1) with high statistics. Superposed, the fit- 



ted form (3.S), as well as the Gaussian and 
exponential components separately. 



Similarly, Fig. Q shows the same method applied to simulations in D = 2. We use a 



2048^ lattice with a^ = 0.4, aX = 0.6. We have 



av 



1.6, d = \plla\i ~ 3.535, and we use 



— = 0.025 — 0.4. Again, the rescaled curves agree very weh, and again (Fig. |5|) the simple 
form ( |3.9| ) is an excellent fit. The fit parameters Oj, /3j are given in Table |l|. 
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Figure 6: The rescaled propagators (| 
£) = 3 for different T. 



Figure 7: One rescaled propagator (r//x = 
0.1) with high statistics. Superposed, the fit- 
ted form ( |3.9| ), as well as the Gaussian and 
exponential components separately. 



Finally, Figs. P and 1^ show results in D = 3, with a fit to the appropriate version of 
dS^gj). We use a 256^ lattice, ap, = 1.0, A = 0.6 to give c?v^ = 10, d = V^. J = 0.08 - 0.64. 
The match is again very good, and we fidn the fit parameters ctj, /3j given in Table |l|. 

For the very large F curve (magenta), the agreement in the UV is not perfect, a result 
of having a large defect density, n ~ 0.08, nd ~ 0.11. This means that 11 percent of all 
lattice points are in the core of a wall, and it is therefore not surprising that our separation 
of scales does not hold. The wall profile is no longer a simple tanh of r jd. But in the IR 
the agreement is still good, and for all smaller damping rates even the UV performs well. 



With the parameters Oj, /3j in (3.10) known, it is possible to determine the defect 
density n from the correlator by a simple fit to ( 3.1C| ), keeping only n as the free parameter. 
We demonstrate this by calculating n as a function of F in this way. Fig. IT] shows the 
"counting" densities n and the "fit" n, resulting from this procedure. The agreement is 
excellent (within 2 percent), except for very small damping rates F in ID, where presumably 
residual thermal noise introduces the discrepancy. Note that for each Z), the calibration was 
performed with a separate, high-statistics, simulation at one particular value of F. These 
are not included in Fig. [l|, and so although because of the successful rescalings Figs. |2[ |^, 
|6| the match could be anticipated, it is a non-trivial result. 



6. Conclusion 



We have shown that the equal-time momentum space two-point field correlator in a defect 
configuration produced by the Kibble mechanism separates in an UV part encoding the 
kink profile (in terms of d) and an IR part encoding the distribution of kinks (in terms of 
the density n) . We carried out numerical simulations to obtain a very good approximation 
to the exact form, which turned out to be the sum of a Gaussian and a position-space 



exponential. This is a somewhat surprising result, since for instance assuming randomly 
distributed kinks in D = 1 yields only the exponential. Furthermore, we found that the 
IR part of the correlator only depends on the defect density and scales with it in a simple 
way. This is a highly non-trivial result and allows us to determine the defect density with 
a one-parameter fit to the field correlator. 

This is important because direct counting of defects, which usually works well in clas- 
sical field theory simulations, is not a meaningful procedure in full quantum field theory 
because the state of the system is not described by a classical field configuration. This can 
be seen concretely in many non-equilibrium quantum field theory techniques, such as the 
2PI formalism, which describes the dynamics in terms of the correlation functions. Our 
results can be employed in such calculations to determine the produced number of defects, 
or to check the validity of these techniques in non-linear situation. Indeed, we used this 



approach in an earlier work |29| to show that at next-to-leading order in the 1/N expansion, 
the 2PI formalism fails to describe defect formation in 1+1 dimensions. 

There is still significant room for improvement in our understanding of the defect 
signature in correlation functions. The present results are only valid at relatively low 
temperature and weak coupling. Thermal and quantum fluctuations change the form of the 
correlator, both by giving a direct contribution to the field correlator and by changing the 
kink profile [^ , and more work is needed to disentangle these effects from the contribution 
due to the kink distribution. Other directions for future work include generalisation to other 
types of defects, such as vortices and monopoles. 
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